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In this Letter we report the first LDA+DMFT (method combining Local Density Approxima- 
tion with Dynamical Mean-Field Theory) results of magnetic and spectral properties calculation 
for paramagnetic phases of FeO at ambient and high pressures (HP). At ambient pressure (AP) 
calculation gave FeO as a Mott insulator with Fe 3d-shell in high-spin state. Calculated spectral 
functions are in a good agreement with experimental PES and IPES data. Experimentally observed 
metal-insulator transition at high pressure is successfully reproduced in calculations. In contrast to 
MnO and Fe2 03 (d'' configuration) where metal-insulator transition is accompanied by high-spin 
to low-spin transition, in FeO (d^ configuration) average value of magnetic moment \/< /i^ > is 
nearly the same in the insulating phase at AP and metallic phase at HP in agreement with X-Ray 
spectroscopy data (Phys. Rev. Lett. 83, 4101 (1999)). The metal-insulator transition is orbital 
selective with only t2g orbitals demonstrating spectral function typical for strongly correlated metal 
(well pronounced Hubbard bands and narrow quasiparticle peak) while Cg states remain insulating. 

PACS numbers: 74.25. Jb, 71.45. Gm 



Introduction- For many years one of the central issues 
of condensed matter physics is the metal-insulator tran- 
sition (MIT) in d- or /-elements compounds [Ij. The 
most spectacular examples are pressure-driven transi- 
tions from wide gap Mott insulators to metallic state for 
transition metal oxides. For MnO and Fe203 {d^ con- 
figuration) metal-insulator transition is accompanied by 
high-spin to low-spin transition (HS-LS). Recently MIT 
in those materials was successfully described theoreti- 
cally by LDA-I-DMFT (method combining Local Density 
Approximation with Dynamical Mean-Field Theory) [2| 
calculations [3|, Ij]. 

Iron oxide also exhibits MIT under high pressure. Re- 
sistivity measurements showed that FeO becomes metal- 
lic at pressures exceeding 72 GPa 5] . Correct description 
of MIT under pressure in wiistite (Fei_a;0) is crucial in 
Earth science because iron oxides are believed to be ma- 
jor constituents of Earth mantle. 

At ambient pressure and room temperature FeO has 
cubic rocksalt Bl structure [6]. Below Neel temperature 
TAr=198 K FeO transforms into rhombohcdral structure 
that could be viewed as a slight elongation along cube di- 
agonal of the original cubic structure. Under pressure at 
room temperature rhombohcdral distortion is observed 
at w 15 GPa and this structure is preserved up to at 
least 140 GPa [TJ, Isj. This transformation to rhombo- 
hcdral structure was believed to accompany long-range 
magnetic ordering due to increasing of Neel temperature 
with pressure [9|. However recent neutron diffraction 
study of wiistite at room temperature under pressure [10| 
showed the absence of magnetic peaks corresponding to 
antifcrromagnetism. At high pressures and temperatures 



(P >120 GPa and T >1000 K) FeO transforms into NiAs 



B8 phase [11 1 



In contrast to MnO and Fe203 it is not clear if FeO 
undergoes HS~LS transition with increase of pressure. 
Controversial experimental evidences were obtained for 
this problem. Mossbauer spectroscopy [l^l shows that 
quadrupole splitting appears between 60 and 90 GPa at 
room temperature. The authors interpreted that as LS 
diamagnetic state. On the other hand high pressure X- 
ray emission spectroscopy |1S] demonstrates that the 
satellite feature in Fe K/3 line associated with HS Fe^'^ 
state does not disappear up to 143 GPa. 

Electronic structure calculations in standard Density 
Functional Theory (DFT) methods predict an antiferro- 
magnetic metallic ground state [lj| in contrast to exper- 
imentally observed insulator with an optical band gap of 
2.4 eV [151] -The LDA-^U method ^Ifi] has been success- 
fully applied to investigate strongly correlated transition 
metal oxides and predicted an insulating ground state 
in FeO at ambient pressure [ITj. Further investigation 
done by Gramsch et al. [18[ for stoichiometric wiistite 
has showed that using the value of Coulomb parameter 
U that reproduces experimentally observed energy gap at 
ambient pressure one can obtain metal-insulator transi- 
tion in LDA+U calculations for unrealistically high pres- 
sures only. 

MIT in transition metal oxides with pressure can 
be successfully described using LDA+DMFT calcula- 
tions [3J, y. In the present work we demonstrate that 
LDA-hDMFT method reproduces MIT for FeO with pres- 
sure. However in contrast to MnO and Fe203 MIT is 
not accompanied by high-spin to low-spin transition and 



metallic spectral function is observed only for t2g orbitals 
while Cg states remain insulating. 

Method.- The LDA+DMFT method [2] calculation 
scheme is constructed in the following way: first, a Hamil- 
tonian Hlda is produced using converged LDA results 
for the system under investigation, then the many-body 
Hamiltonian is set up, and finally the corresponding self- 
consistent DMFT equations are solved. The calculations 
presented below have been done for crystal volumes cor- 
responding to values of pressure up to 140 GPa and room 
temperature. Since no structure transition has been ob- 
served at low temperatures [8] and NiAs phase appears 
above 1000 K only all calculation were performed for sim- 
ple NaCl (Bl) cubic crystal structure with lattice con- 
stant scaled to give a volume corresponding to applied 
pressure [5|. Ab-initio calculations of electronic struc- 
ture were obtained within the pseudopotential plane- 
wave method PWSCF, as implemented in the Quan- 
tum ESPRESSO package [l9J. Hamiltonians Hlda in 
Wannier function (WF) basis 20|, |21| were produced us- 
ing projection procedure that is described in details in 
Ref.li. 

The WFs are defined by the choice of Bloch functions 
Hilbert space and by a set of trial localized orbitals that 
will be projected on these Bloch functions. The basis 
set includes all bands that are formed by 0-2p and Fe- 
3d states and correspondingly full set of 0-2p, and Fe-3d 
atomic orbitals to be projected on Bloch functions for 
these bands. That would correspond to the extended 
model where in addition to d-orbitals all p-orbitals are 
included too. 

The resulting 8x8 p — d Hamiltonian to be solved by 
DMFT has the form 
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where U^^ 

the occupation number operator for the d electrons with 
orbitals a or /3 and spin indices a or ct' on the i-th site. 
The term Hdc stands for the d-d interaction already ac- 
counted for in LDA, so called double-counting correction. 
In the present calculation the double-counting was chosen 
in the following form Hdc = ^^('^■dmft — ^)i- Here ndmft is 
the self-consistent total number of d electrons obtained 
within the LDA-I-DMFT, U is the average Coulomb pa- 
rameter for the d shell and / is unit operator. 

The elements of U^o matrix are parameterized by U 
and Jh according to procedure described in [23|. The 
values of Coulomb repulsion parameter U and Hund ex- 
change parameter Jh were calculated by the constrained 
LDA method 24 1 on Wannier functions [2^. Obtained 
values Jh=0-89 eV, U = 5 eV are close to previous 
estimations [l8[. The effective impurity problem for 
the DMFT was solved by the hybridization expansion 
Continuous-Time Quantum Monte-Carlo method (CT- 
QMC) [25|. Calculations for all volumes were performed 




n ' r 

-2 2 

Energy [eV] 



FIG. 1: (Color online) Spectral function of Fe-d states vs. 
pressure obtained in LDA-fDMFT (CT-QMC) calculations 
at room temperature. 



in the paramagnetic state at the inverse temperature 
P = 1/T = 40 eV-i corresponding to 290 K. Spectral 
functions on real energies were calculated by Maximum 
Entropy Method (MEM)!!!]. 

Results and discussion- The Fe d band is split by crys- 
tal field in triply degenerated t2g and doubly degenerated 
eg subbands. LDA fails to describe insulating ground 
state of FeO at AP and for all volumes FeO is metallic. 

Including Coulomb correlation effects in frames of 
LDA+DMFT method results in high spin state wide gap 
Mott insulator for ambient pressure phase (APP) of FeO 
in agreement with experimental data. The calculated 
energy gap value of about 2 eV agrees well with IPES 
measurement 27| value 2.5 eV and optical spectrum [15| 
value 2.4 eV. The occupation numbers for Fe d orbitals 
are n(e(,)=0.54 and n{t2g)~0-68. The average value of lo- 



cal magnetic moment •\/< fil > is B.S/ib. Those numbers 
agree very well with high-spin state of Fe+^ ion (d^ con- 
figuration) in cubic crystal field: 2 electrons in eg states 
(n(eg)=l/2) and 4 electrons in t2g states (n(i2g)=2/3) 
with magnetic moment value 4 fj.B- Spectral functions 
A{u!) for all pressure values calculated by MEM using 
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FIG. 2: (Color online) Total spectral function of FeO in am- 
bient pressure phase calculated within LDA+DMFT (CT- 
QMC) (/3=40 eV~^) (solid blue line) in comparison with com- 
bined PES and IPES experimental data (red dots) from 

Ref. Ei,!!!. 



Green function G{t) from CT-QMC calculations are pre- 
sented in the Fig. [1] The spectral function for ambient 
pressure phase (APP) shows well defined insulating be- 
havior for all d-orbitals. However the energy gap for eg 
states is nearly two times larger than for t2g states indi- 
cating that the latter orbitals are closer to MIT than the 
former ones. Figure [5] contains calculated total spectral 
function compared with_spectrum combined from PES 
and IPES experiments 



M 291. The theoretical and ex- 



perimental curves are in a good agreement. 

LDA-I-DMFT calculation made for small volume values 
corresponding to high pressures gave metallic state for 
FeO (see Fig.JJI) starting from 60 GPa in agreement with 
experiment [5|]. One can see that t2g orbitals become 
metallic whereas Cg ones remain insulating. This behav- 
ior reminds the orbital selective Mott transition (OSMT) 
in ruthenates [30|. Occupation number values in Fe-d 
shell are practically not changed comparing with APP 
and are n(eg)=0.55 n(i2g)=0.68 at 140 GPa. The mag- 
netic moment value decreases on a few percent only and 
is 3.5/xb at 140 GPa. The only interpretation for those 
values is that an iron d-shell in high pressure metallic 
phase of FeO still corresponds to high-spin state of Fe"^^ 
ion. This conclusion agrees well with analysis of high 
pressure X-Ray emission spectroscopy experiment made 
in Ref. [13[. The occupation numbers and magnetic mo- 
ment vs. pressure are presented in the Fig. [31 One can 
see that all curves exhibit the kink at 60 GPa. We argue 
that this feature is due to MIT and corresponding re- 
construction of spectral function at Fermi level. Spectral 
functions A{uj) for t2g in the Fig. [1] for pressure values 
larger then 60 GPa become typical for strongly correlated 
metal close to MIT: well pronounced Hubbard bands and 
narrow quasiparticle peak. A(cu) for Cg is still insulat- 
ing with Hubbard bands only but energy gap value is 



FIG. 3: (Color online) Magnetic moments (black squares) and 
occupancies of t2g (red cicles) and eg (blue triangles) shells vs. 
pressure obtained in LDA+DMFT (CT-QMC) calculations. 



strongly decreased comparing with APP (see Fig.[T]). 

To understand these results the following simple model 
was used. The model has two semicircle DOS of the same 
width with 3 orbital and 4 electrons. One orbital is non- 
degenerate and two other orbitals are degenerate. The 
centers of gravity and DOS widths were taken from ab 
initio LDA calculations. In this model non-degenerate 
orbital stands for e^-orbital in FeO and two others for 
t2g. Occupations in model in HS state are 1/2 for non- 
degenerate orbital (the same as in realistic LDA-I-DMFT 
calculation for FeO) and 3/4 for degenerate orbitals com- 
paring with 2/3 in the case of t2g orbitals in FeO. The 
Kanamori parametrisation of Coulomb repulsion (with 
the same U=5 eV and J=0.89 eV) was used. Note, that 
corresponding matrix elemetns U^'Z (eq. [l} are set to 
be the same for all orbitals. The model was solved using 
DMFT (CT-QMC) method and obtained spectral func- 
tions for two values of pressure (APP and 140 GPA) are 
presented in the Fig. |4l 

The orbital selective metal-insulated transition 
(OSMT) was reproduced in these calculations. Since 
DOSes for all three orbitals have the same width (in 
contrast to OSMT [30(] in ruthenates where two bands 
have very different widths) and actual structure of DOS 
is neglected we can conclude that effects of different 
degeneracy of orbitals and deviation from half filling are 
the driving force of this separate transition. It is known 
that critical value of Coulomb interaction parameter 
Uc needed for metal-insulator transition in half-filled 
degenerate Hubbard model is Uc ~ vNU^^^ — NJ 
[3l| (iV is degeneracy and U^^^ is critical U value 
for non-degenerate case). That means that for more 
degenerate i2g orbitals one needs larger effective U 
value to become insulating than for less degenerate eg 
orbitals. In addition to that for half-filled states an 
estimation for effective Ueff value is C/-I- {N — 1) J while 
for the occupancy one electron more then half-filling 
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FIG. 4; (Color online) Left panels - LDA densities of 
states (DOS) (solid lines) and correponding model semicir- 
cle DOS(daslied lines). Right panels - spectral functions from 
model DMFT (CT-QMC) calculations for two values of pres- 
sure. Non degenerate orbital (black lines) reproduces eg- 
orbitals and 2 times degenerate one reproduces t2g orbitals 
(red line). 



an estimation is f/e// = U — J. Then for 2/3 filled t2g 
orbitals one needs much larger U value to drive them 
into insulating state than for half-filled eg states. 

In conclusion- We have performed LDA-I-DMFT cal- 
culation for FeO at room temperature and values of pres- 
sure from the ambient one till 140 GPa. In the agree- 
ment with experiment spectral function for FeO at AP 
demonstrates an energy gap of about 2 eV. At the pres- 
sures higher then 60 GPa FeO is metallic but only for 
t2g orbitals while Cg states remain insulating that corre- 
sponds to orbital selective Mott transition scenario. The 
MIT obtained in our calculations is not accompanied by 
change of spin state and FeO has HS with large local mo- 
ment in APP and all HPP. This result agrees with high 
pressure X-Ray emission spectroscopy data. 
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